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In the context of f{R) theories of gravity, we address the problem of finding static and spherically 
symmetric black hole solutions. Several aspects of constant curvature solutions with and without 
electric charge are discussed. We also study the general case (without imposing constant curvature). 
Following a perturbative approach around the Einstein-Hilbert action, it is found that only solutions 
of the Schwarzschild-(Anti) de Sitter type are present up to second order in perturbations. Explicit 
expressions for the effective cosmological constant are obtained in terms of the f{R) function. Fi- 
nally, we have considered the thermodynamics of black holes in Anti-de Sitter space-time and found 
that this kind of solutions can only exist provided the theory satisfies Ro + f{Ro) < 0. Interestingly, 
^\ I this expression is related to the condition which guarantees the positivity of the effective Newton's 

, constant in this type of theories. In addition, it also ensures that the thermodynamical properties 

' in f(R) gravities are qualitatively similar to those of standard General Relativity. 
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I. INTRODUCTION 



In the last years, increasing attention has been paid to modified theories ofgravity in order to understand several 
open cosmological questions such as the accelerated expansion of the universe and the dark matter origin • Some 
O of those theories modify General Relativity by adding higher powers of the scalar curvature R, the Riemann and 
Ricci tensors or their derivatives Q. Lovelock and f{R) theories are some examples of these attempts. It is therefore 
5-H , quite natural to ask about black holes (BH) features in those gravitational theories since, on the one hand, some BH 
I signatures may be peculiar to Einstein's gravity and others may be robust features of all generally covariant theories 

of gravity. On the other hand, the results obtained may lead to rule out some models which will be in desagreement 
CNJ ■ with expected physical results. For thoses purposes, research on thermodynamical quantities of BH is of particular 
J> I interest. 

. In this work we will restrict ourselves to the so called f{R) gravity theories (see 'A]) in metric formalism in Jordan's 
' frame. In this frame, the gravitational Lagrangian is given by f{R) where f{R) is an arbitrary function of R and 
, Einstein's equations are usually fourth order in the metric (see Q for several proposed /(i?) functions compatible with 
local gravity tests and other cosmological constraints) . An alternative approach would be to use the Einstein's frame, 
where ordinary Einstein's gravity coupled to a scalar plus a massive spin-2 field is recovered. Even if a mathematical 
correspondence could be established between those two frames, in the last years some controversy has remained about 
their physical equivalence. 

• • ' Previous literature on /(i?) theories fd\ proved in Einstein's frame that Schwarzschild solution is the only static 
. 5^ ' spherically symmetric solution for an action of the form R + aR^ in Z? = 4. In [7] uniqueness theorems of spherically 
, symmetric solutions for general polynomial actions in arbitrary dimensions using Einstein's frame were proposed (see 
^ • also [1| for additional results). See also [§| for spherical solutions with sources. 

5^ I Using the euclidean action method (see for instance [l^, [IH) in order to determine different thermodynamical 
quantities, Anti de Sitter (AdS) BH in /(i?) models have been studied p^. In [l^ the entropy of Schwarzschild-de 
Sitter BH was calculated for some particular cosmologically viable models in vacuum and their cosmological stability 
was discussed. 

BH properties have been also widely studied in other modified gravity theories. For instance, [3 studied BH 
in Einstein's theory with a Gauss-Bonnet term and cosmological constant. Different results were found depending on 
the dimension D and the sign of the constant horizon curvature k. For k — 0,-1, the Gauss-Bonnet term does not 
modify AdS BH thermodynamics at all (only the horizon position is modified with respect to the Einstein-Hilbert 
(EH) theory) and BH are not only locally thermodynamically stable but also globally preferred. Nevertheless for 
fc = +1 and D = 5 (for D > 6 thermodynamics is again essentially that for AdS BH) there exist some features not 
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present in the absence of Gauss-Bonnet term. Gauss-Bonnet and/or Riemann squared interaction terms were studied 
in [iGj concluding that in this case phase transitions may occur with k — —1 . 

Another approach is given by Lovelock gravities, which are free of ghosts and where field equations contain no 
more than second derivatives of the metric. These theories were studied in 17] and the corresponding entropy was 
evaluated. 

The paper is organized as follows: in section 2 we present some general results for f{R) gravities for interesing 
physical situations in metric formalism. In sections 3 and 4, BH in /(i?) gravities are studied and explicit Einstein's 
field equations are presented for static and spherically symmetric metrics. Section 5 is devoted to find perturbative 
solutions for static and spherically symmetric background metric: general metric coefficients are found depending on 
/(i?) derivatives evaluated at background scalar curvature. Sections 6 and 7 are widely devoted to study thermody- 
namical quantities and their consequences in local and global stability for some particular f{R) models. Finally, we 
include some conclusions. 



II. GENERAL RESULTS 



In order to study the basics of the solutions of general f(R) gravity theories, let us start from the action 

S ^ Sg + S^ (1) 

where Sg is the D dimensional gravitational action: 

= f d^'^VWl (R + fiR)) (2) 



IGttGi _ 

with Gjj = being the D dimensional Newton's constant, Md the corresponding Planck mass, g the determinant 

of the metric qab, {A, B — 0,1, D~ 1), R the scalar curvature and R+f{R) is the function defining the theory under 
consideration. As the simplest example, the EH action with cosmological constant A^i is given by f{R) = —{D—2)A]j. 
The matter action Sm defines the energy-momentum tensor as: 

J.AB ^ 2 

V I .9 I ^9AB 

From the above action, the equations of motion in the metric formalism are just: 

RAsil + fiR)) - \iR + f{R)) QAB + (VaVb - gABn)f{R) + SttGdTab = (4) 

where Rab is as usual the Ricci tensor and □ — V^V"^ with V the usual covariant derivative. Thus for the vacuum 
EH action with cosmological constant we have: 

1 D — 2 

Rab ~ -^RgAB H ^ — ^ogAB = (5) 

which means Rab = Ad9ab and R = DAjj. Coming back to the general case, the required condition to get constant 
scalar curvature solutions R = Rq (from now Rq will denote a constant curvature value) in vacuum implies: 

RABil + .nR))-l9AB{R + f{R))^0 (6) 

Taking the trace in previous equation, Rq must be a root of the equation: 

2(1 + /'(i?o)) Ro~D {Ro + f{Ro)) = (7) 

For this kind of solution an effective cosmological constant may be defined as A^-^ = Ro/D. Thus any constant 
curvature solution R — Rq with 1 -I- /'{Rq) ^ fulfills: 

Rq + f{Ro) 
^■^^-2(l + /'(i?o))^^^ 

On the other hand one can consider: 

2R (1 + fiR)) -DiR + fiR)) = (9) 
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as a differential equation for the /(i?) function so that the corresponding solution would admit any curvature R value. 
The solution of this differential equation is just: 

f{R) = ai?^/2 - R (10) 

where a is an arbitrary constant. Thus the gravitational Lagrangian becomes proportional to aR^/^ which will have 
solutions of constant curvature for arbitrary R. The reason is that this action is scale invariant since a/Go is a 
non-dimensional constant. 

Now we will address the issue of finding some general criteria to relate solutions of the EH action with solutions 
of more general ,f{R) gravities, not necessarily of constant curvature R. Let qab a solution of EH gravity with 
cosmological constant, i.e.: 

1 D — 2 

Rab - ^RgAB + -^Ad9ab + SttGdTab = (11) 
Then qab is also a solution of any f{R) gravity, provided the following compatibility equation 

f'{R)RAB - IgAB [/(i?) + {D^ 2)Ad] + (Va Vb - gABa)f'{R) = (12) 

obtained from ^ is fulfilled. In the following we will consider some particularly interesting cases. The simplest 
possibility is obviously vacuum (Tab = 0) with vanishing cosmological constant Ajj = 0. Then the above equation 
PTjl becomes: 

Rab = ^RgAB (13) 

which implies R = and Rab — 0. Consequently gAB is also a solution of any f{R) gravity provided /(O) — 0, 
which is for instance the case when f{R) is analytical around R = 0. When the cosmological constant is different 
from zero (A^ 7^ 0), but still Tab = 0, we have also constant curvature with Rq — DAo and Rab — AoQab- 
Then the compatibility equation (fT2|) reduces to ([7]). In other words, qab is also a solution of f{R) provided 
f{DAo) = Ad{2 — D + 2f'{DAD))- Notice that it would also be a solution for any i?o in the particular case 
f{R) = ai?^/2 - R. 

Next we can consider the case with A^ = and conformal matter (T = T^ = 0). For a perfect fluid this means 
having the equation of state p = {D — l)p where p is the pressure and p the energy density. In this case (jll[) implies 

R = ; Rab = SttGdTab (14) 

Then, provided /(O) — /'(O) — 0, gAB is also a solution of any f{R) gravity. This result could have particular 
interest in cosmological calculations for ultrarelativistic matter (i.e. conformal) dominated universes. For the case of 
conformal matter with non vanishing Ad we have again constant R = Rq with Rq — DAjj and gAB is a solution of 
f{R) provided that once again /{DAd) = Ad{2 - D + 2f'{DAD)). 



III. BLACK HOLES IN f{R) GRAVITIES 

Now we consider the external metric for the gravitational field produced by a non rotating object in /(i?) grav ity 
theories. The most general static and spherically symmetric D > A dimensional metric can be written as (see |18|): 

ds2 = e-2*('-)A(r)d<2 - A-\r)dr^ - r^dn]j_^ (15) 

or alternatively 

ds2 = A(r)d<2 - p-^{r)dr^ - r^<^^D-2 (16) 

where dri|)_2 is the metric on the S^^^ sphere and identification A(r) e^^*'^'")yl(r) and p{r) = A{r) can be 
straightforwardly established. 

For obvious reasons the <I>(r) function is called the anomalous redshift. Notice that a photon emitted at r with 
proper frequency ujq is measured at infinity with frequency lOoq = e~*'^'')-y/yl(r)wo- As the metric is static, the scalar 
curvature Rm D dimensions depends only on r and it is given, for the metric parametrization (fT51) . by: 

R{r) = \[D^ -^D + Q + rA'{r){-2D + ir<^'{r)+4) 

- r'^A"{r) - A{r) {D^ - 5D + 2r^'i>'{rf - 2{D - 2)r$'(r) - 2r^'i>"{r) + 6)]. (17) 
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where the prime denotes derivative with respect to r. At this stage it is interesting to ask about which are the most 
general static and spherically symmetric metrics with constant scalar curvature i?o- This curvature can be found 
solving the equation R = Rq. Then it is inmediate to see that for a constant = (f>o the general solution is: 

A{r) = 1 + a,r^-^ + a,r'-^ D^^f ^^^^ 

with ai and 02 being arbitrary integration constants. In fact, for the particular case D — R^^ — Q and $0 ~ 0, the 
metric can be written exclusively in terms of the function: 

+ V + (19) 

By establishing the identifications ai = —2GnM and 02 = Q^, this solution corresponds to a Reissner-Nordstrom 
solution, ie. a charged massive BH solution with mass M and charge Q. Further comments about this result will be 
made below. 



IV. CONSTANT CURVATURE BLACK-HOLE SOLUTIONS 



By inserting the metric (|15p into the general /(i?) gravitational action 5*^ in and making variations with respect 
to the A{r) and <i>(r) functions, we find the equations of motion: 



(2 -D)il + f'iRWir) - r [r{R)R'{rf + f"{R){<^' {r)R! [r) + i?"(r))] = 



(20) 



and 



2rA{r)f"'{R)R'{rf + f"{R)[2DA{r)R'{r) - AA{r)R'{r) + 2rA{r)R"{r) + A'{r)rR'{r)] + 
g' {R)[-2rA{r)<i>' {rf + 2DA{r)^'{r) - 4A(r)$'(r) - rA"{r) + 2rA{r)^"{r) + 

A'{r){2 -D + 3r$'(r))] - r{R + f{R)) = (21) 
where /', /" and /"' denote derivatives of ,f{R) with respect to the curvature R. 

The above equations look in principle quite difficult to solve. For this reason we will firstly consider the case of 
constant scalar curvature R = Ro solutions. Then the equations of motion reduce to: 



(2 - i^) (1 + /'(i?))$'(r) = 



and 

R+f{R) + {l + f'{R)) 



A"{r) + {D- 2)^^ - {2D - 4)^1M£M _ 3A'(r)$'(r) + 2A(r)$'2(r) - 2A(r)$"(r) 



As commented in the previous sections, the constant curvature solutions of /(i?) gravities are given by: 

Df{R^) 
" 2{l + nRa))-D 

whenever 2(1 + /'(i?o)) ^ -D. Thus from ^ $'(r) = and then ([23]) becomes 

A'{r) 



Ro + f{Ro) + (1 + f{Ro)) 

Coming back to ([25|). and using ([24|). we get 



A"{r) + {D-2y- 



= 



A"{r) + iD^2)^ = -^Ro 
r D 



(22) 

= 
(23) 

(24) 
(25) 
(26) 



This is a /(i?)-independent linear second order inhomogeneous differential equation which can be easily integrated to 
give the general solution: 



A[r) = Ci + Car' 



,3-r> 



D{D-l) 



(27) 
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which depends on two arbitrary constants Ci and C2. However this solution has no constant curvature in the general 
case since, as we found above, the constant curvature requirement demands Ci ~ I. Then, for negative Rq, this 
solution is basically the D dimensional generalization obtained by Witten [ll| of the BH in AdS space-time solution 
considered by Hawking and Page With the natural choice $0 = the solution can be written as: 

^M = l--|-r + 77- (28) 



where 



with 



= ^^^^ 

being the area of the D — 2 sphere, P = —D{D — l)/i?o is the asymptotic AdS space scale squared and M is the mass 
parameter usually found in the literature. 

Thus we have concluded that the only static and spherically symmetric vacuum solutions with constant (negative) 
curvature of any /(i?) gravity is just the Hawking-Page BH in AdS space. However this kind of solution is not the 
most general static and spherically symmetric metric with constant curvature as can be seen by comparison with the 
solutions found in (|18p . Therefore we have to conclude that there are constant curvature BH solutions that cannot be 
obtained as vaccum solutions of any f{R) theory. As we show below, in the D = A case, we see that the most general 
case can be described as a charged BH solution in /(i?)-Maxwell theory. 

Indeed, let us consider now the case of charged black holes in f{R) theories. We will limit ourselves to the D = A case, 
since in other dimensions the curvature is not necessarily constant. The action of the theory is now the generalization 
of the Einstein-Maxwell action: 

Sg = T^ [ d''xV\J]{R + f{R)-F,,,Fn (31) 



where F^t^ = dfj,A^ — d^A^. Considering an electromagnetic potential of the form: = {V{r),0) and the static 
spherically symmetric metric (|15|) . we find that the solution with constant curvature Rq reads: 

Vir) = ^ 

AM = = 1 - ^ + (i±mO)Q! _ Ro^. (32) 

r 12 

Notice that unlike the EH case, the contribution of the black-hole charge to the metric tensor is corrected by a 
(1 + /'(i?o)) factor. 



V. PERTURBATIVE RESULTS 



In the previous section we have considered static spherically symmetric solutions with constant curvature. In EH 
theory this would provide the most general solution with spherical symmetry. However, it is not guaranteed this to be 
the case also in f{R) theories. The problem of finding the general static spherically symmetric solution in arbitrary 
f{R) theories without imposing the constant curvature condition is in principle too complicated. For that reason 
in this section we will present a perturbative analysis of the problem, assuming that the modified action is a small 
perturbation around EH theory. 

Therefore let us consider a f{R) function of the form 

fiR) =-{D-2)AD + ag{R) (33) 

where a ^ 1 is a dimensionless parameter and g{R) is assumed to be analytic in a. By using the metric parametrization 
given by (jl6p the equations of motion become: 

A(r)(l + f'iR)){2^l{r)[iD-2)X'ir)+rX"ir)]+rX'ir)^/{r)} 

- 2X{rf {2^,{r)[{D - 2)R'{r)f"{R) + rf^^\R)R'{rf + rR"{r)r{R)] + rR' {r)^i' {r)f" {R)] 

- rM(r)A'(r)2(l + /'(i?)) + 2rX{rf{R + f{R)) = (34) 
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- A(r)At'(r) [2(1? - 2)A(r) + rA'(r)] (1 + /'(i?)) 

+ ^l{r)[2X{r)R'{r) [2(0 - 2)A(r) + rA'(r)] /"(i?) + r(l + /'(i?))(A'(r)2 - 2A(r)A"(r))} 

- 2rX{rf{R + f{R)) = (35) 

where prime denotes derivative with respect to the corresponding argument and R = R{r) is given by ()17p . Now, 
assuming that the A(r) and /i(r) functions appearing in the metric (I16p are also analytical in a, they can be written 
as follows 



A(r) = Ao(r)+^a'A,(r) 

oo 

fi{r) = ^o(f) + 



(36) 



i=l 



where {Ao(r), //o('')} are the unperturbed solutions for the EH action with cosmological constant given by 



Aio(0 = 1 + 



Ci Ad 2 



,D-3 



(D-l) 



Ao(0 = -C2(i?-2)(i?-l)Aio(r) 



(37) 



which are the standard BH solutions in a I? dimensional AdS spacetime. Note that the factor C2 can be chosen by 
performing a coordinate t reparametrization so that both functions could be indentified. For the moment, we will keep 
the background solutions as given in ([37]) and we will discuss the possibility of getting A(r) — fi{r) in the perturbative 
expansion later on. 

By inserting ([33]) and p6|) in (|34p and ((35)) we obtain the following first order equations: 



(D - 3)Aii(r) + r^i(r) + = ^ 



(38) 



Ci(2?-3)r3-^ + -^r2 



C2 [Ci(D - l)r'^-^ - Az3r2 + i? - l] g{R,y + 
+ C72(i?-2)(i?-l)(Az5r2-i? + 3)/ii(r) 

[2^2(1 -i?)r2 Ads' (i?o)+rA'i(r)] - 



Ai(0 



1 + Cir' 



whose solutions are: 



Ai(r) = C^{B-\){p-2) 



{CxC^-G^G^){p~2){D-\) 



rD-3 



[C^D - 2)Ad + C2 ig{Ro) - 2ADg'iRo))] : 



f . C3 ^ (g(i?o) - 2Ang'iRo)) , 



rD-3 {D-2){D-l) 



Up to second order in a the equations are: 



(n , , (.g(flo) - 2Acff'(j?o)) ( 2D . ^ 2 

[D - 3)^2 (r) + r^2(^) H -R — 9 [Ro] ~ -j^ — T^^og [Ro) r 



D -2 



D -2 



(39) 



(40) 
(41) 

(42) 
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D - 1 



A2(r) + C2{D - 2){D - 1) (-Ai^r^ + - 3) /i2(r) 



4-D 



D 



jj X'^ir) - ^3^4(1? -2)iD- 1) (-Az5r2 + - 3) 



„3--D 



C2 [(i^-l)(Cir3-^ + l)-Az3r2] 



2A,,g'(i?o)'+g(i?o) 



2, ,^ , /2i?A,3.g"(i?o) „„A 4DA|3g'(7?o)g"(i?o) 



5'(i?o) 



C4[Ci{D-l] 



„3--D 



2][2Az,5'(i?o)-5(i?o)]r' = 



whose solutions are: 



A2(r) = C6 + 



CeCi + (C3C4 - C2C5)(i5 - 2){D - 1) 



+ {giRoj - ^^og [Ro)) + 623 (i?o) - ■ 



£> - 1 



D-2 



D-2 



(43) 



(44) 



/"2(0 = 



C5 , (5(i?o) - 2ADg'{Ro)) {2DADg"{Ro) - {D - 2)g'{Ro)) 



,,D-3 



{D-2)^{D-1) 



(45) 



Further orders in a^-"^'--- can be obtained by inserting previous results in the order 3,4,... ones to get 
{A3,4,...(r),/L(3^4_...(r)} but of course the corresponding equations become increasingly complicated. 

Notice that from the obtained results up to second order in a, the corresponding metric has constant scalar 
curvature for any value of the parameters Ci, C2, . . . , Cg. As a matter of fact, this metric is nothing but the standard 
Schwarzschild-Ad5' geometry, and can be easily rewritten in the usual form by making a trivial time reparametrization 
as follows: 



A(r) = X{r) [-C2{D^ + 3D ~ 2) + C4 {D^ - 3D + 2) a + Cea^ + 0{a^)] 
-p{r) = n{r) 



(46) 



Therefore, at least up to second order, the only static, spherically symmetric solutions which are analytical in a 
are the standard Schwarzschild-AdS' space-times. 

On the other hand, taking the inverse point of view, if we assume the solutions to be of the AdS BH type at any 
order in the a expansion we can write: 



A(r) EE ^(r) = 1+ f +Jr^ 
as solution for the Einstein equations ([34]) and ((35|) with the gravitational lagrangian ([33| and 



(47) 



Rs = Rs + S^iCiO; 



SOO 7^4 



{D-l) 

where Rs and Ci are arbitrary constants and the Ji coefficients can be determined from ([7]): 

R-{D- 2)Kd + ag{R) + 2{D - 1) J(l + ag' {R)) = 



(48) 



(49) 



with R — ~D{D — 1) J. Expanding previous equation in powers of a it is possible to find a recurrence equation for 
the Ji coefficients, namely for the J; (with I > 0) coefficient, we find: 



i-i 



{2-D){D-l)Ji + Y, E 



1 



i— cond.l 
l-l l-k-1 



2iD-l)Y^J, E E 



fc=0 i=0 



'id. 2 



n!i2! . . .i(-fc-i! 



{Jiy'{j2Y' ■ ■ ■ {Ji-k-iy'-''-'g^'+^\Ro) 



(50) 
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with Rq = —D{D — 1) Jo = DAd, where the first sum is done under the condition 1 given by: 

/-I i-i 

im — i, hn e NU {0} and ^ mim — I — ^ (51) 



m— 1 m — 1 

and the second one under the condition 2: 

i-fc-l l-k-l 



m—1 m—1 



For instance we have: 



im = i, im G N U {0} and ^ mim = I — k — 1 (52) 
^(.g ; Az,) 



Ji 



(i?-2)(i^-l) 



, _ A(g; i^, Az,)[(i^ - 2)ff^(Eo) - 2DADg"{Ro)] 

"^^ ~ (iP-2)2(i^-l) ^^^> 

where ^(5 ; Z?, A^) = g{Ro) - 2ADg'{Ro). 

Now we can consider the possibihty of removing Ajj from the action from the very beginning and stiU getting an 
AdS BH solution with an effective cosmological constant depending on g{R) and its derivatives evaluated at Rq = 0. 
In this case the results, order by order in a up to order a^, are: 

MAn^O) = 



{D-2){D-1) 

As we see, in the context of /(i?) gravities, it is possible to have a BH in an AdS asymptotic space even if the initial 
cosmological constant A^i vanishes. 

To end these two sections, we can summarize by saying that in the context of /(i?) gravities the only spherically 
symmetric and static solutions of negative constant curvature are the standard BH in AdS space. The same result 
applies in the general case (without impossing constant curvature) in perturbation theory up to second order. However, 
the possibility of having static and spherically symmetric solutions with non constant curvature cannot be excluded 
in the case of f{R) functions which are not analytical in a. 



VI. BLACK-HOLE THERMODYNAMICS 



In order to consider the different thermodynamic quantities for the f{R) black- holes in AdS, we start from the 
temperature. In principle there are two different ways of introducing this quantity for the kind of solutions we are 
considering here. Firstly we can use the definition coming from Euclidean quantum gravity p^ . In this case one 
introduces the Euclidean time t — it and the Euclidean metric ds'^ is defined as: 

- ds| = -d<7^ - r^dnl,_2 (55) 

where: 

da^ = e-2*MA(r)dT2 + A''\r)dr^. (56) 

The metric corresponds only to the region r > ru where th is the outer horizon position with A{rH) = 0. Expanding 
d(T^ near th we have: 

da^ ^ e-2*('^-)A'(rH)pdr2 + (57) 

A'{rH)p 

where p — r — ru- Now we introduce the new coordinates R and 9 defined as: 

9 - ie-*('--)A'(rff)r 
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so that: 

da^ = R^AO^ +dR^. (59) 

According to the Euchdcan quantum gravity prescription r belongs to the interval defined by and (3e = ^/Te- On 
the other hand, in order to avoid conical singularities, 6 must run between and 27r. Thus it is found that 

Te = -^e-*(^-U'(r^) (60) 
47r 

Another possible definition of temperature was firstly proposed in [20| stating that temperature can be given in 
terms of the the horizon gravity /C as : 

T. . ^ (61) 

where K. is given by: 

IC = hm ^i^. (62) 

'■^^^ V\9ttgrr\ 

Then it is straightforward to find: 

Tk^Te. (63) 

Therefore both definitions give the same result for this kind of solution. Notice also that in any case the temperature 
depends only on the behaviour of the metric near the horizon but it is independent from the gravitational action. By 
this we mean that different actions having the same solutions have also the same temperature. This is not the case 
for other thermodynamic quantities as we will see later. Taking into account the results in previous sections and for 
simplicity we will concentrate only on constant curvature AdS BH solutions with $ = as a natural choice and: 

f>D-3 2 



Then, both definitions of temperature lead to: 



/^-Vr= .^ (65) 



Notice that the temperature is a function of rn only, i.e. it depends only on the BH size. In the limit th going to 
zero the temperature diverges as T ~ and for rn going to infinite T grows linearly with rH- Consequently T 

has a minimum at: 



^HO=l\/^^-^ (66) 



corresponding to a temperature: 



The existence of this minimum was established in [Tq] for Z? = 4 by Hawking and Page long time ago and it is well 
known. More recently Witten extended this result to higher dimensions [ll|. The minimun is important in order to 
set the regions with different thermodynamic behaviors and stability properties. For Z3 = 4, an exact solution can be 
found for rn- 

2/3 



2'/' 9^ + V 12 + 81^ _ (24)1/3 



TH = I ' , 1/3 (68) 

62/3 [ QiM. + Ji2 + 8i:f 
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Thus, in the Rs <C I hmit, we find rn = Rs, whereas in the opposite case I ^ Rs, we get rn = {PRsY^"^- For the 
particular case D = 5, rjj can also be exactly found to be: 

A p2 \ 

(69) 

which goes to i?| for Rs <C I and to IRg for I ^ Rg. Notice that for any T > Tq, we have two possible BH sizes: 
one corresponding to the small BH phase with rn < rno and the other corresponding to the large BH phase with 

In order to compute the remaining thermodynamic quantities, the Euclidean action 

1 f .D^ 



SO that: 



^i^ = -Y^^y d^2;Vgi(i? + /(i?)) (70) 

is considered. When the previous expression is evaluated on some metric with a periodic Euclidean time with period 
/?, it equals f3 times the free energy F associated to this metric. Extending to the f{R) theories, the computation 
by Hawking and Page , generalized to higher dimensions by Witten [ll[ , we compute the difference of this action 
evaluated on the BH and the AdS metric which can be written as: 

AS, = , ^o + /(^") Ay (71) 
where Rq = —D{D — and AF is the volume difference between both solutions which is given by: 

AV=|^(/Vr^-.r^) (72) 

Notice that from this expression it is straightforward to obtain the free energy F. We see that provided — (i?o + 
f{Ro)) > 0, which is the usual case in EH gravity, we have > for rjj < ^ and F < for rn > I- The temperature 
corresponding to the horizon radius rn — I will be denoted Ti and it is given by: 

Notice that for D > 2 we have Tq <Ti. 

On the other hand, the total thermodynamical energy may now be obtained as: 

j,_dASE_ [Rq + fiRo))MP 

^ - ~W " W^) ^ ' 

where M is the mass defined in This is one of the possible definitions for the BH energy for f{R) theories, 

see for instance [2l| for a more general discussion. For the EH action we have f{R) = — (D — 2)Ku and then it is 
immediate to find E = M . However this is not the case for general f{R) actions. Notice, that positive energy in AdS 
space-time requires i?o + f{Ro) < 0. Now the entropy 5" can be obtained from the well-known relation: 

S = I3E- PF. (76) 

Then one gets: 

(Ro + f{Ro))PAD^2{rH) 
^ 8{D-1)Gd 

where Ao~2{'''h) is the horizon area given by Ao-2i'''H) = f^^^ IJ.d-2- Notice that once again positive entropy 
requires i?o -I- /(-Rq) < 0. For the EH action we have i?o + f{Ro) = — 2(1) — and then we get the famous 

Hawking-Bekenstein result [2^ 



AD-2irH) 

^Gd 



(78) 
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Finally we can compute the heat capacity C which can be written as: 

dE drn 



dE 



drn dT 



Then it is easy to find 



C 



-{Ro + f{Ro)){D - 2)MD-2rg"^P jP - lyj, + {D- W 

iD^l)rj,-iD-3)P 



For the particular case of the EH action we find: 

(iJ - 2)f,n^,r^-^ (D - \)r\ + (D - 3)P 
4Gd {D-l)rl-(D-i)p- 

In the Schwarzschild limit I going to infinity this formula gives: 



C 



AG 



< 



D 



(79) 



(80) 



(81) 



(82) 



which is the negative well-known result for standard BH. In the general case, assuming like in the EH case (i?o + 
f{Ro)) < 0, we find C > for vh > tho (the large BH region) and C < for r// < rno (the small BH region). For 
fH ~ tho {T close to To) C is divergent. Notice that in EH gravity, C < necessarily implies F > since Tq < 7i. 

In any case, for f{R) theories with i?o + /(-Ro) < 0> we have found an scenario similar to the one described in full 
detail by Hawking and Page in [lo| long time ago for the EH case. 

For T < Tq, the only possible state of thermal equilibrium in an AdS space is pure radiation with negative free 
energy and there is no stable BH solutions. For T > Tq we have two possible BH solutions; the small (and light) 
BH and the large (heavy) BH. The small one has negative heat capacity and positive free energy as the standard 
Schwarzschild BH. Therefore it is unstable under Hawking radiation decay. For the large BH we have two possibilities; 
ii Tq < T < Ti then both, the heat capacity and the free energy are positive and the BH will decay by tunneling into 
radiation, but if T > Ti then the heat capacity is still positive but the free energy becomes negative. In this case the 
free energy of the heavy BH will be less than that of pure radiation. Then pure radiation will tend to tunnel or to 
collapse to the BH configuration in equilibrium with thermal radiation. 

In general f{R) theories one could also in principle consider the possibility of having Rq + /(i?o) > 0. However 
in this case the mass and the entropy would be negative and therefore in such theories the AdS BH solutions would 
be unphysical. Therefore Rq + /(-Ro) < can be regarded as a necessary condition for f{R) theories in order to 
support AdS BH solutions. Using (Ojthis condition implies 1 + /'(-Ro) > 0. This last condition has a clear physical 
interpretation in f{R) gravities (see [23[ and references therein). Indeed, it can be interpreted as the condition for 
the effective Newton's constant Geff = G_d/(1 + f'{Ro)) to be positive. It can also be interpreted from the quantum 
point of view as the condition which prevents the graviton from becoming a ghost. 



VII. PARTICULAR EXAMPLES 



In this section we will consider some particular /(-R) models in order to calculate the heat capacity G and the free 
energy F as the relevant thermodynamical quantities for local and global stability of BH's. For these particular models, 
Rq can be calculated exactly by using ([7]). For the sake of simplicity we will fix the /^-dimensional Schwarzschild 



radius in 



as R 



D-3 



— 2. The models we have considered are: 



A. Model I: f{R) = a{-Rf 
Substituting in ([7]) for arbitrary dimension we get 

We will only consider non-vanishing curvature solutions, thus we find: 



-Ro — — 



2-D 
{2(3 -D)a 



l/(/3-l) 



(83) 



(84) 
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Since D is assumed to be larger than 2, the condition (2/3 — D)a < provides weU defined scalar curvatures Rq. Two 
separated regions have thus to be studied: Region 1 {a < 0, /3 > D/2} and Region 2 {a > 0, (3 < D/2}. For this 
model we also get 

Notice that in Region 1,1 + /'{Rq) > for D > 2, since in this case /3 > 1 is straightforwardly accomplished. In 
Region 2, we find that for D > 2, the requirement Rq + f{Ro) < 0, i.e. 1 + /'(i?o) > 0, fixes (3 < 1, since this is the 
most stringent constraint over the parameter (3 in this region. Therefore the physical space of parameters in Region 
2 is restricted to he {a > 0, (3 < 1}. 

In Figs. 1-3 we plot the physical regions in the parameter space {a, (3) corresponding to the different signs of (C, F). 



B. Model II: f{R) = -(-i?)" exp(g/J?) - R 
In this case, a vanishing curvature solution appears provided a > 1. In addition, we also have: 

To get i?o < the condition q{2a ~ D) < must hold and two separated regions will be studied: Region 1 {q > 
0, a < D/2} and Region 2 {q < 0, a > D/2}. 

In Figs. 4-6 we plot the regions in the parameter space {a, q) corresponding to the different signs of (C, F). 



C. Model III: f{R) = R{logaRy - R 

A vanishing curvature solution also appears in this model. The non trivial one is given by 

Ro = ^exp (^) (87) 

Since Rq has to be negative, a must be negative as well, accomplishing aRo > and since aR, and therefore aRo, has 
to be bigger than one to have a positive number powered to q, what imposes g > as can be read from the argument 
of the exponential in the previous equation. Therefore there exists a unique accessible region for parameters in this 
model: a < and q > 0. 

In Figs. 7-8 we plot the regions in the parameter space {a, q) corresponding to the different signs of (C, F). 



D. Model IV: f{R) = ~a^^^ 



l+/3(f )" 

This model has been proposed in [24] | as cosmologically viable. Throughout this section, we consider n ~ 1 for this 
model. Hence imposing /'(i?o) = e we get 

. - -^^^^ m 

hence a relation between ci , D and e can be imposed and therefore this model would only depend on two parameters 
a and p. A vanishing curvature solution also appears in this model and two non trivial curvature solutions are given 
by: 



± a [(ci ~2)D + A± ^VciD'' - 8D + 16^ 

= TT^^TR i^^) 



2/3(1} -2) 

The corresponding 1 + f'{Ro) values for ([55)1 are 



1 + /'(i?±) = 1 - , , ^ (90) 



(VciL)^ -8L» + 16 ± y/^D) 
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where ci > and ci > {8D — 16)/D^ are required for real Rq solutfons. Since 1 + /'(i?o) > is required, that means 
that sign(i?^) = sign(a/3). It can be shown that 1 + f'{R^) is not positive for any allowed value of ci and therefore 
this curvature solution Rq is excluded for our study. 

1 + /'{Rq) > only requires ci > for dimension D > 4 and therefore e < is required according to (|88p . 
Therefore only two accesible regions need to be studied: Region 1 {a > 0, /3 < 0} and Region 2, {a < 0, /3 > 0}. 

In Figs. 9-10 we plot the thermodynamical regions in the parameter space (a,/3) for a chosen e = — 10~®. Note 
that 1 + /'{Rq) does depend neither on a nor on /3 and that Rq only depend on the quotient a/f3 for a fixed ci. 

VIII. CONCLUSIONS 

In this work we have considered static spherically symmetric solutions in f{R) theories of gravity in arbitrary 
dimensions. After discussing the constant curvature case (including charged black-holes in 4 dimensions), we have 
studied the general case without imposing, a priori, the condition of constant curvature. We have performed a 
perturbative analysis around the EH case which makes possible to study those solutions which are regular in the 
perturbative parameter a. We have found explicit expressions up to second order for the metric coefficients, which 
give rise to constant curvature (Schwarzschild AdS) solutions as in the EH case. 

On the other hand, we have also calculated thermodynamical quantities for the AdS black holes and considered 
the issue of the stability of this kind of solutions. We have found that the condition for a f{R) theory of gravity to 
support this kind of black holes is given by Rq + f{Ro) < where Rq is the constant curvature of the AdS space-time. 
This condition has been seen to imply also that the effective Newton's constant is positive and that the graviton does 
not become a ghost. For these f{R) gravities the qualitative thermodynamic behavior of the BH is the same as the 
one found by Hawking and Page for the AdS BH but the value of some thermodynamic magnitudes is different for 
different f{R) gravities. 

Finally we have considered several explicit examples of f{R) functions and studied the parameter regions in which 
BH in such theories are locally stable and globally preferred, finding the same qualitative behaviour as in standard 
EH gravity. 
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(a) Model I, D = 4, Region 1, a < 0, /3 > 2. 
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(b) Model I, D = i, Region 2, a > 0, /3 < 1. 
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(b) Model I, D = b, Region 2, o > 0, /3 < 1. 



Figure 2: Thermodynamical regions in the (a,/3) plane for Model I in D = 5. Region l(left), Region 2 (right). 
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(b) Model I, D = 10, Region 2, a > 0, /3 < 1. 



Figure 3: Thermodynamical regions in the (a,/3) plane for Model I in D = 10. Region l(left). Region 2 (right). 
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(a) Model II, D = 4, Region 1, a < 2, 5 > 0, 
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(b) Model II, D = 4, Region 2, a > 2, g < 0. 



Figure 4: Thermodynamical regions in the {a,q) plane for Model II in D = 4. Region l(left), Region 2 (right). 
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Model II D=5, Region 1 
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(a) Model II, D = 5, Region 1, « < 2.5, g > 0. 
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(b) Model II, D = 5, Region 2, a > 2.5, q < 0. 



Figure 5: Thermodynamical regions in the (a, g) plane for Model II in D = 5. Region l(left), Region 2 (right). 
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(b) Model II, D = 10, Region 2, a > 5, g < 0. 



Figure 6: Thermodynamical regions in the (a, g) plane for Model II in D = 10. Region l(left). Region 2 (right). 
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(a) Model III, D = 4, a < 0, q > 0. 



(b) Model III, D = 5, a < 0, g > 0. 



Figure 7: Thermodynamical regions in the {a, q) plane for Model III in D = 4 (left) and D — 5 (right). 
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Figure 8: Thermodynamical regions in the (a, q) plane for Model III in D = 10. 
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Figure 9: Thermodynamical regions in the {\a\, \l3\) plane for Model IV in D = 4 (left) and D — 5 (right). 
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Figure 10: Thermodynamical regions in the (|a|, \/3\) plane for Model IV in Z) = 10. 



